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1 Introduction 



Our knowledge of the 7-ray sky has increased dramatically during this last decade, due principally 
to the 7-ray instruments on board the Compton Gamma Ray Observatory (CGRO)[|l[|: EGRET, 
a spark chamber plus calorimeter instrument with 7-ray sensitivity in the energy interval 30 
MeV to 30 GeV; COMPTEL, a Compton telescope in the interval 0.1 to 30 MeV; BATSE, an 
omnidirectional x-ray and soft 7-ray "burst" detector consisting of large Nal scintillators sensitive 
to 30 keV to 2 MeV photons (with smaller spectroscopic Nal crystals for measurements up to 110 
MeV); and OSSE, consisting of Nai-Csl phoswiches detecting photons of 0.1 to 10 MeV. Their 
observations have revolutionized our understanding of such extragalactic phenomena as blazars 
and gamma ray bursts (GRBs), as well those within our own Galaxy such as pulsars. 

A salient characteristic of these extragalactic sources is their extreme time variability. In 
the case of blazars, flare states are observed to occur in which the luminosities vary by over 
an order of magnitude within time-scales as small as hours. GRBs are singular, extragalactic 
events due either to the the coalescence of massive compact objects or to the hypernovae of 
massive early stars, whose luminosities decay over periods of seconds to days depending upon the 
wavelength observed. For such objects, the importance of correlated, multiwavelength observations 
is paramount. By monitoring the relative time lags of blazar flare peaks at radio, optical, x-ray, 
and 7-ray wavelengths, for example, the properties of the source regions at the bases of AGN 
(active galactic nucleus) jets are constrained. Gamma ray observations of blazars by EGRET have 
been particularly crucial, as the luminosity of these objects peaks in the GeV-to-TeV part of the 
electromagnetic spectrum. The transient nature of these sources therefore argues for continuous 
coverage of the 7-ray sky. 

Until recently, observations with CGRO and other space-based and ground-based telescopes 
have provided coverage of these and other sources up to the limiting sensitive energy of EGRET, 
approximately 30 GeV. From there, a gap in our knowledge of the 7-ray sky spectrum has ex- 
isted up to 200-300 GeV, where ground-based 7-ray shower detectors presently have their energy 
thresholds. It is possible that within this gap there are novel features in the 7-ray sky, such as 
a gamma-ray line or continuum emission from postulated neutralino annihilation at the center 
of the Galaxy. Future instruments with sensivity in this unexplored region, such as AGILE ||], 
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recently approved by the Italian Space Agency (ASI), or the proposed GLAST M payload, may 
uncover exciting new phenomena. 

At the present time, however, our view of the 7-ray sky has diminished. With the effective 
turnoff of the EGRET 7-ray instrument on the CGRO due to the nearly complete consumption 
of its spark chamber gas, there will soon be no operating instrument capable of observing the 
7-ray sky in the energy interval ~ to ~ 10 2 GeV. Ground-based 7-ray detectors, based on 
the atmospheric Cerenkov technique (ACT)Q, turn on at current energy thresholds of ~200-500 
GeV. Within the next several years energy thresholds for some ACT observatories are expected to 
go to as low as 20 to 50 GeV, but lower thresholds than these are unlikely to be achieved due to the 
sizable effect of Earth's magnetic field on the 7-ray-induced air showers, and the lower Cerenkov 
photon yield which must be detected against the night sky background. For 7-ray energies much 
lower than 10 2 GeV, then, spaced-based detectors are required. This observational gap for the 
energy window ~ 10 _1 to 10 2 GeV will continue to exist for the next several years. Eventually, 
this gap will be eliminated with the launch of a next-generation 7-ray satellite, such as GLAST 
P], but such a mission is unlikely to occur before the year 2005. Also the AGILE satellite ||, to 
be launched in 2002, will have a limited sensitivity above 50 GeV. 

In this paper we describe how the Alpha Magnetic Spectrometer (AMS) can largely fill this 
gap by acting as a 7-ray detector with sensitivity in the energy interval of 0.3 to 100 GeV during 
its three-year mission on board the International Space Station Alpha (ISSA) from 2003 to 2006. 
AMS, described in detail elsewhere, has as its primary mission the search for cosmic ray antinuclei, 
which are detected via the negative curvature of their trajectories through a magnetic spectrometer 
(see Sect. ||). 

It is of interest, however, to analyze if a large particle detector like AMS, operated on the Space 
Station for three years starting in 2003, could be used to detect high energy 7-rays . There are 
at least two ways this could be done without significantly affecting the experiment sensitivity to 
antimatter: a) adding a light (e.g. 0.3Ao) converter at the entrance of the magnetic spectrometer, 
either passive (e.g., a high-Z thin plate) or active (e.g. a multilayer tracking detector), or b) 
implementing an high granularity imaging shower detector at the bottom of the experiment. 
Option (a) consists in a minor modification to the baseline instrument, while option (b), although 
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would give a much higher 7-ray detection capabilities, requires the addition of a new detector. In 
this paper we present a study which has been perfomed on option (a). We will show that in this 
option, AMS can also detect 7-rays with performance characteristics similar to EGRET in the 
energy region of 0.3 to 20 GeV, and with significantly enhanced capabilities between 20 and ~ 100 
GeV, a region which is not well unexplored. We refer to this modified instrument as "AMS/7". We 
show that AMS/7 can continue the valuable work of EGRET by providing continued monitoring 
of extragalactic and galactic 7-ray sources and by participating in multiwavelength observational 
campaigns. In addition, AMS/7 wu l have unprecedented sensitivity to the 7-ray sky between 
the energies of 20 to ~ 100 GeV (albeit at a level somewhat lower than required for detection 
of known sources with power law spectra), so that AMS/7 might provide us with unexpected 
discoveries in this region. 

The next section gives a brief description of the baseline AMS instrument, and the modi- 
fications required for adding 7-ray detection capability in option (a). Section 3 describes the 
performance characteristics of AMS/7 determined from Monte Carlo analyses, and Section 4 uses 
these results to address the question of what 7-ray astrophysics can be done with AMS. The 
results obtained analyzing option (a) give a hint of the physics potential which could be reached 
implementing option (b) for the flight of AMS on the Space Station; option (b) would gain a 
factor of about 5 in 7-rays statistics over option (a). Detailed analysis of option (b), however, wil 
be the subject of future work. 

2 The Alpha Magnetic Spectrometer 
2.1 The Baseline Instrument 

The Alpha Magnetic Spectrometer has been built by a large international collaboration of high 
energy physics institutions from the U.S., Italy, China, Finland, France, Germany, Taiwan, Russia, 
and Switzerland. It recently had a successful test flight on the Space Shuttle mission STS-91 in 
June, 1998 where it was carried in the cargo bay and observed for several days in both the 
zenith and nadir directions, the latter for measurement of albedo cosmic ray backgrounds. It 
is scheduled to be secured to an external payload attachment point on the International Space 
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Table 1: AMS silicon tracker parameters (in parenthesis the values used on the precursor flight). 



Number of planes 


6 


Accuracy (bending plane) 


10 fim 


Accuracy (non bending plane) 


30 fxm 


Number of channels 


163936 




(58368) 


Power consumption 


400 W 




(180 W) 


Weight 


130 kg 


Silicon Area (double sided) 


5.4 m 2 




(2.4 m 2 ) 



Station Alpha in 2003, where it will remain as a zenith-pointing instrument for a minimum of 3 
years. 

The basic design of AMS is shown in Figure |l|. The magnet spectrometer consists of a per- 
manent ring dipole magnet made of very high grade Nd-B-Fe rare earth material whose magnetic 
energy product and residual induction are respectively (.£?i^) max > 50 x 10 6 G-Oe and 14,500 G, 
yielding an analyzing power of (BL 2 ) = 0.14 T-m 2 with less than 2 tonnes of magnet mass. 

Four silicon strip detector tracking planes are located within the magnetic volume, with a 
fifth and sixth plane located just above and below the magnet. These are based on technology 
developed for the Aleph and L3 vertex detectors at LEP @ , and consist of 300/im thick, double- 
sided, high purity silicon wafers. The total area of this tracker is about an order of magnitude 
larger than the size of previous microstrip silicon detectors. The momentum resolution of the 
silicon spectrometer is given in Figure || at low rigidities (below 8 GV) the resolution is dominated 
by multiple scattering (Ap/p ~ 7%), while the maximum detectable rigidity (Ap/p ~ 100%) is 
about 500 GV. The parameters of the silicon spectrometer are given in Table 1. In addition to 
measuring particle rigidity, the silicon planes will provide six independent measurements of dE/dx 
for charge determination. 
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Four time-of-flight (ToF) scintillator planes (two above and two below the magnet volume) 
measure the particle velocity with a resolution of 120 ps over a distance of 1.4 m. The ToF 
scintillators also measure dE/dx, allowing a multiple determination of the absolute value of the 
particle charge. 

A solid state Cerenkov detector below the magnet provides an independent velocity measure- 
ment. In addition, a scintillator anticounter system is located within the inner magnet wall, 
extending to the ToF scintillators. 

The performance of this instrument as a nuclear cosmic ray detector, and its sensitivity to an 
antinucleus cosmic ray component, is discussed elsewhere ||. 

We should also note that by the time AMS is attached to the International Space Station it 
may have undergone significant changes from the baseline design considered here. In particular, 
the permanent magnet may be replaced by a superconducting magnet, which would consider- 
ably improve AMS's 7-ray detection performance (although it would increase the 7-ray energy 



threshold as well; see 3.3). Since the detector is still undergoing significant design changes from 
the baseline instrument flown on the Space Shuttle, we have chosen to fix on that design which 
currently exists as integrated hardware. The addition of other components, such as a transition 
radiation detector, a segmented calorimeter, a solid Cerenkov radiator, and/or a superconducting 
magnet will only improve the performance of AMS/7 described in this paper. 

2.2 Modifications of AMS for 7-ray Detection 

In option (a), the conversion of AMS to AMS/7 requires the addition of two hardware components: 
(1) A high-Z converter plate (e.g. tungsten) of 0.3 radiation lengths (x = 0.3Ao), which converts 
7-rays into electron-positron pairs. The momenta of each electron and positron are measured 
in the spectrometer, and their weighted sums give the primary photon's energy and incidence 
direction. (2) A stellar attitude sensor gives the angular orientation of AMS with respect to the 
celestial sphere to an accuracy of better than 1.5 x 10~ 4 radians. The stellar attitude sensor is 
required because the Space Station cannot provide attitude information to better than ~ 2°, while 
the AMS/7 angular resolution approaches 1' for E"f > 100 GeV. Commercial units, such as Ball 
Aerospace Systems Division's CT-633 Stellar Attitude Sensor, provide angular resolutions to 0.5' 
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with very modest weight and power consumption (2.5 kg and < 10 W). 

Determination of the converter thickness and position is based on the following considerations: 

(1) The 7-ray converter must not degrade the sensitivity of AMS for antimatter detection. There- 
fore the converter mass cannot be distributed uniformly throughout the tracking volume, which 
would otherwise provide optimal angular and energy resolution. Instead, the converter must lie 
as a single plate just above the second TOF scintillator, S2. 

(2) The loss of nuclear antimatter candidates to nuclear interactions in the converter plate must 
be insignificant: At x = 0.3Xq, this amounts to only ~ 2% of cosmic-ray 4 He being lost. 

(3) The sensitivity for detection of astrophysical point sources of gamma rays is nearly indepen- 
dent of x for x <C Xq for background-limited (as opposed to statistics-limited) detection. This 
sensitivity degrades as x/Xq — ► 1. 

(4) The aperture for rare events and for measurement of the diffuse gamma ray background in- 
creases with x. 

(5) The fractional energy loss of the electron pair due to bremsstrahlung within the converter is 
oc x/Xq, degrading energy resolution as x increases. 

(6) Can the presence of the converter plate create background that will degrade the cosmic ray 
He/He target sensitivity of 10 -9 ? By requiring clean tracks in the spectrometer volume we do not 
expect converter-generated contamination, as the converter lies above the magnet volume. This 
is corroborated both by the SHuttle flight data (where an He/He of ~ 10~ 6 was obtained |7T| ) 
and by Monte Carlo analyses. This question is also being examined using helium and carbon data 
from GSI at Darmstadt. 

For this anaysis we choose the value x = 0.3Xo, for which point source sensitivity is still 
optimal, the energy resolution is acceptable, and nuclear interaction losses are negligible. 

3 The Performance of AMS/7 
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3.1 Monte Carlo Calculations 

A full-instrument GEANT Monte Carlo code was run to determine the performance of AMS/7. 
Gamma rays with fixed energies, ranging from 0.3 to 100 GeV, were thrown isotropically at the 
detector over an opening angle of 50°. 

For those 7-rays which converted in the detector (most within the converter plate) the charged 
particles were tracked through the inhomogeneous field of the ring dipole permanent magnet. The 
inhomogeneous field used was that measured for the AMS magnet that was recently flown on the 
Space Shuttle (flight STS-91), with a field strength increased by 7% to match the expected field 
strength of a second magnet which might be constructed for the Space Station mission. 

All the physical processes for electrons and 7-rays were "on" in the GEANT code. Bremsstrahlung 
photons of energies < 20 MeV were not followed, although all bremsstrahlung energy losses 
were included. The hit positions in the silicon tracker planes were assigned gaussian errors with 
a = 15/zm and 30^m for the bending and nonbending planes, respectively. (Test beam mea- 
surements give a bending plane Si tracker resolution below 10 /im 0.) Reconstruction of the 
particle trajectories were performed with a quintic spline fit routine ||] to determine their mo- 
menta and charge sign. Assignment of a hit position within the 6-plane Si tracker to either the 
converted electron or positron was done for all possible combinations (64 in all) of assignments; 
that combination whose trajectory reconstruction x 2 was minimum was selected as the event. All 
secondary charged particles were tracked in the detector volume until they stopped or exited from 
the detector mother volume. 

The 7-ray candidate events were those that passed the hardware trigger of SI • S2- S3 -54, where 
51 to 5*4 are the ToF counters (the converter plate lies between SI and S2). If reconstruction 
of particle trajectories occurred successfully, further cuts were placed on these candidate events 
before acceptance as a 7-ray event: 

1) the reconstructed 7-ray zenith angle # 7 must be < 50°; 

2) the reconstructed 7-ray energy E 1 must be > 0.15 GeV; 

3) there are at least two reconstructed tracks; 

4) the reconstructed energy of the most energetic particle be > 0.15 GeV; 

5) the number of anticoincidence panels hit (there are 16 such panels along the interior wall of 
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the magnet) be < 4. 

Primary 7-ray energy and incidence direction were determined by adding the fitted momenta 
vectors of all secondaries evaluated at the converter plate to obtain the primary momentum 
vector. (More sophisticated weighting methods were found to have negligible effect on 7-ray 
angular resolution.) The results of the Monte Carlo runs are discussed in the subsections below. 

3.2 Point Source Sensitivity 

Astrophysical sources of 7-rays fall into two categories: point sources, such as blazars and GRBs, 
and diffuse sources, such as the cosmic isotropic 7-ray background and possibly 7-ray emission 
from neutralino annihilation near the center of our Galaxy. Some galactic sources fall into an 
intermediate category of "extended" sources, such as pulsar nebula and supernova remnants, 
where the source angular size may be larger than the point spread function of the 7-ray detector. 

For point sources, a detector's key figure of merit is its point source sensitivity, defined as the 
minimum source flux required to achieve a specified level of detection significance. The significance 
S of a detection is given schematically by 

N(> Et) 



(i) 



where N(> E t ) and B(> E t ) are respectively the total number of detected photons from the 
source and the number of background photons falling within the source area above an energy E t . 
E t is usually, but not always, the instrument's threshold energy. N and B are a function of the 
effective detection area A{E) of the instrument, its angular resolution expressed as a solid angle 
£l(E), the viewing time t, and the differential source and background spectra dN/dE and 0IB/0IE: 



r°° dN 

N(>E t ) = J —A(E)tdE, (2) 



and 

ro ° dB 



B(> Et) = r ^-A(E)Q(E)t dE. (3) 
J Ex dhj 

All extragalactic objects (blazars) seen by EGRET are well-represented by power law spectra 



dN ( E V 7 
dE =U0 



llGevJ (- 2 -GeVr\ (4) 
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where 7 is the differential spectral index of the source, and no is the differential flux at 1 GeV. We 
define our point source sensitivity as that value of no which gives a 5a signal integrated over the 
period of one year of operation of AMS. It is evident that the point source sensitivity is a function 
of the source's spectral index; for these calculations we have assumed a source differential spectral 
index of 7 = 2, which is close to the mean of the source spectral indices of blazars observed by 
EGRET §. 

The background flux dB/dE contains components from the isotropic extragalactic 7-ray back- 
ground radiation and from the galactic diffuse radiation (the latter due to the decay of 7r°s pro- 
duced in collisions of cosmic rays with the interstellar medium). The former has recently been 



determined by EGRET to be [10 



jd / p \2.10±0.03 

(7.32 ± 0.34) X 10- (g^) (c^.s.sr.GeV)-. ,5) 

The galactic diffuse radiation is a strong function of position in the sky, and been measured and 
modelled by the EGRET group [11]. We have used here their data to generate background spectra 
j[g(E, 5, a) over 3° bins in declination 5 and right ascension a. 

The solid angle Q(E) over which the background must be integrated when viewing a source is 
given by Q(E) = iraQ 8 (E), where o^s, following EGRET's usage, is defined as the angular radius 
within which 68% of the source photons fall. The converter thickness of 0.3 Xq largely dominates 
all multiple scattering effects, so that the angular and energy resolution of reconstructed primary 
photons will be completely dominated by multiple Coulomb scattering (MCS) and bremsstrahlung 
energy losses of the electrons within the converter plate. Figure |3] confirms that the Monte 
Carlo result for AMS/7's angular resolution is very close to that estimated from MCS within the 
converter plate: 

/ TP \ -0.956 

a 6 A 8 MS (£) = 0.88° (-^) , (6) 



which is to be compared to EGRET's angular resolution [12] of 

/ p \ -0.534 
= 1.71° (^JL.) . (7) 

The effective area A{E) of the instrument is also a function of the zenith angle 6 of the incident 
photon, as shown in Figure [3 for 3 GeV photons. Although the dependence of the conversion 
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probability cancels that of the projected area, the requirement that S3 and £4 ToF scintillators 
be hit causes the effective area to drop to zero for > 50°. An integration of effective area 
A(E, 0) over solid angle gives instrument aperture, shown in Figure |5l Below a 7-ray energy of 
0.3 GeV, the converted electrons begin to have too small a radius of curvature to escape from 
the magnet volume, and detection efficiency plummets. Above 100 GeV the converted electron 
and positron often do not spatially diverge beyond the two-hit resolution distance of the silicon 
trackers, causing significant deterioration in 7-ray energy resolution. These then define the limits 
of the energy window for AMS/7. 

Because AMS will be rigidly attached to the ISSA which itself will be in 51.6° inclination orbit, 
AMS will not spend equal amounts of time viewing all directions in the celestial sphere. However, 
AMS/7 has finite detection area out to 9 = 50°, so that AMS/7 does view the entire sky. Full 
sky coverage is obtained roughly every three months due to the precession of the orbital plane 
of the Space Station about Earth's pole [13|. Because of this precession, the viewing area x time 
product over the sky is only a function of E and declination 5, and not of right ascension a. By 
following the orbit of AMS over one year, at each time element determining the AAt increment 
for each position in the sky viewable by AMS at that time, a total At(E, 5) areaxtime product is 
obtained, and is shown in Figure |(| 

With the above results and Eqs. AMS/7's 5a point source sensitivity can be obtained. 



The significance is calculated using the likelihood ratio of Li and Ma [ 14 1 , 



where a = t on /t g, t on being the time spent viewing the source, and t Q g being the time spent 
viewing off-source, and N on and N q q are the total integrated on-source (signal plus background) 
and off-source (background) photon counts. Since at 1 GeV roughly 10~ 3 of the instrument's 
field-of-view (FOV) is contained in the <768 error circle, we take a = 10 -3 . For any given position 
in the sky, that computed value of no which gives S = 5 at 1 GeV is the point source sensitivity 
for that sky position. Figure [?] shows the one-year point source sensitivity for a lower cutoff energy 
Et of 1 GeV. 

One cannot directly compare the point source sensitivity of AMS/7 to that of EGRET, since 
AMS is not a pointable instrument. EGRET achieves a flux sensitivity of I m in(> 0.1 GeV) « 
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10 cm - 2 s with a 2-week viewing period. Since it took EGRET one year to map the full sky, 
where each sky segment was viewed for roughly 2 weeks, we can compare EGRET's sensitivity 
with that achieved by AMS/7 after one year of operation. By assuming a source differential 
spectrum of E~ 2 , we can convert from EGRET's definition of sensitivity (in terms of integral flux 
above 0.1 GeV) to that of AMS/7 (a differential flux above 1 GeV). In our units the EGRET 5a 
flux sensitivity is 1 x 10 -8 . Over most of the sky, AMS/7's mean sensitivity (no) is estimated to 
be about a factor of 2 lower than that of EGRET. 

3.3 Energy Resolution 

The 7-ray energy resolution is dominated by bremsstrahlung losses in the converter plate, which 
have an average value of l^f"} ~ \ (gx^cose) ~ 0-16/ cos 9 for a converter thickness r = 0.3AV 
In the reconstructed 7-ray energy distributions for primary energies of 0.3, 3, 20, and 100 GeV, 
shown in Figure ||, one sees the effect of bremsstrahlung losses in the converter plate in the large 
low-energy tails of the distributions. The nearly identical appearances of these distributions is a 
consequence of the fact that the probability for the electron to radiate a fraction v of its energy 
as a bremsstrahlung photon is solely a function of v. There are two noticeable deviations from 
this universality in the figure: The cutoff of events below 0.1 GeV is due to the gyroradii of 
the electrons becoming smaller than the magnet field volume dimensions; the broadening of the 
distribution at 100 GeV is due to the fact that trajectory reconstruction errors are beginning to 
dominate bremsstrahlung losses in the energy determination. Figure |9| gives the mean fractional 
energy loss and the rms energy resolution as a function of primary 7-ray energy, and shows the 
deviations from constant values below 1 GeV and near 100 GeV. 

That the reconstructed energy distributions have a nearly "universal" form has important 
consequences for the determination of blazar energy spectra. If a source energy spectrum follows 
a power-law energy distribution, E~^ , and if the energy resolution function of the detector is of the 
form f(E', E) = E~ 1 g(E' /E), where E is the true incident energy, E' is the reconstructed energy, 
and g(x) is an arbitrary function of x, then the reconstructed source spectrum is also a power-law 
spectrum with the same spectral index. There will therefore be little bias in our reconstructed 
spectral indices. This is demonstrated in Section 4.1.2] (Figure 11) where the measured spectra of 
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several blazars, a convolution of the primary 7-ray energy spectrum E 7 and the energy resolution 
function f(E',E), are seen to match the true spectra. 

3.4 Comparison of AMS/7 with EGRET 

The table below summarizes the performance characteristics of AMS as a 7-ray detector. In 
particular, by its comparison to EGRET, one sees that the two instruments perform similarly in 
many respects, one major difference being the energy windows: AMS's energy window is shifted 
up by roughly one order of magnitude from that of EGRET, thereby providing an improved view 
of the sky in the region E 1 ~ 100 GeV. 





AMS 


EGRET 


technique 


magnetic spectrometer 


spark chamber plus calorimeter 


energy window (GeV) 


0.3 to 100. 


0.03 to 30. 


peak effective area (cm 2 ) 


1300 


1500 


angular resolution 


0.77°(E/lGeV)-°- 96 


l.71°(E/lGeV)- - 53A 


half-area zenith angle 


~ 30° 


~ 20° 


total viewing time (yr) 


~3 


~2 


attitude capability 


fixed 


movable 


flux sensitivity 
(ph/cm 2 -s-GeV at 1 GeV) 


~ 0.5 x 1(T 8 


~ 1.0 x 1(T 8 



A few of the table entries require some amplification: (1) The 30 GeV upper limit listed for 
EGRET is due in part to the effects of electromagnetic backsplash from their calorimeter into 
their hardwired anticoincidence scintillator shield that surrounds the instrument. Above ~GeV 
energies this self-induced veto reduces EGRET's aperture roughly as E~°- 5 . Although the EGRET 
collaboration quotes 30 GeV as their maximum energy, in fact they have small but finite aperture 
up to and beyond 100 GeV. (2) The total viewing time for EGRET is listed as 2 years. This takes 
into account the fact that, on average, during approximately one-third of the viewing time the 
pointing direction is occulted by Earth. (3) Unlike EGRET, AMS/7 is n °t a pointing instrument; 
its vertical axis is always aligned with the zenith. Therefore AMS/7 cannot point to targets of 
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opportunity. 

Note that the point source sensitivity is nearly the same for both EGRET and AMS, with 
that of AMS being somewhat (factor ~2) lower (given the caveat of Section |3.2| ), This is because 
EGRET and AMS have similar effective detection areas and angular resolutions. The fact that 
AMS's 7-ray energy threshold is an order of magnitude larger than EGRET's, thereby implying 
an integrated point source flux about an order of magnitude lower for AMS than for EGRET, is 
compensated for by the fact that AMS's angular aperture (half-angle area of 30°) is larger than 
EGRET's (~ 20°), and by the fact that AMS spends 100% of its time pointing to the sky (being 
attached to a gravity-gradient stabilized Space Station), while EGRET typically points one-third 
of the time to Earth. AMS lacks the low-energy end of EGRET's range due to the curvature 
of the electron-positron pair in AMS's magnetic field, which limits the 7-ray aperture to > 300 
MeV. On the other hand, EGRET lacks the high-energy end of AMS's range due to the effect of 
electromagnetic backsplash in EGRET's Nal calorimeter which vetos most events above 30 GeV. 

The above table shows that AMS/7 and EGRET have very similar 7-ray detection capabilities, 
so that AMS will be able to continue and extend the investigation of galactic and extragalactic 7- 
ray sources initiated by EGRET. One main difference between the two detectors, that of significant 
aperture above 30 GeV for AMS/7, may lead to the observation of new phenomena in this 
relatively uncharted region. A more detailed discussion of the astrophysics goals of AMS/7, an d 
in particular the prospect for novel discoveries, is given in the next section. 

4 Gamma Ray Astrophysics Goals of AMS/7 
4.1 Blazar Astrophysics 

Until the launch of CGRO, only one extragalactic source of > 0.1 GeV 7-rays was known, the 
blazar 3C 273, detected by the early 7-ray satellite COS-B EGRET has subsequently in- 

creased the number of detected blazars to over 50 M, and in fact all of the extragalactic sources 
seen by EGRET are blazars. The term blazar is used to describe a wide variety of objects, such 
as flat-spectrum radio quasars, optically violent variable (OVV) quasars, and BL Lac objects, 
all of which are characterized by high luminosity and rapid time variability. Within the "unified 
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model" |T(| these 7-ray-loud sources have a common source of energy: the accretion of matter 
onto a supermassive (1O 7_1O M0) black hole, which in turn funnels much of its expelled energy 
along the black hole spin axis, forming relativistic "jets" of outflowing plasma |17| . When the axis 
of such a jet points towards Earth, relativistic beaming effects dramatically increase the apparent 
luminosity and time variability of these sources, producing a "blazar" signature [^] 

Nearly all EGRET blazars are found to have their luminosity peaking in the 7-ray region, i.e., 
at levels higher than the radio, IR, optical, or x-ray regions flT9|] . Blazars detected by EGRET 
are found to belong to the blazar subclasses of flat-spectrum-radio-quasars and radio-selected 
BL Lacs H, ^, whose 7-ray emissions extend to the tens of GeV, but apparently fall off at 



higher energies, since they are not detected by ground-based instruments whose energy thresholds 
are as low as 200 GeV. Ground-based instruments, on the other hand, have detected a total of 5 
blazars worldwide, all of one subclass, x-ray-selected BL Lacs. Of these 5 blazars, EGRET has 
detected only one (Mrk 421). Clearly an interesting transition is occurring between the effective 
EGRET cutoff of ~ 20 GeV and the present ground-based detector threshold of ~ 200 GeV, 
the observation of which may help determine the 7-ray production mechanisms within the blazar 
jets. It is not known whether the spectral cutoffs at ~ 10 2 GeV of EGRET blazars are due to 
intrinsic processes within the jets, or are due to extinction of the higher energy 7-rays off the 
diffuse, extragalactic soft photon (IR-optical) background j2^] (via 77 — > e + e~ interactions) as 
they propagate over cosmic distances to Earth. 

4.1.1 Source Counts 

The number of blazars that we expect to detect with AMS/7 can be estimated from the EGRET 
blazar population. During its first six years of operation (1991 to 1997) EGRET has claimed the 
detection of 51 blazars |J. Different identification methods, however, have resulted in an EGRET 
blazar count ranging from 42 pi to roughly 60 pi. Of the 71 unidentified 7-ray sources detected 



by EGRET, 14 are at galactic latitudes \b\ > 30° p3j, making them likely blazar candidates. 



Figure 10 shows the flux distribution of these sources (the identified blazars, plus the 14 high- 
latitude unidentified sources). In this figure the flux is the integral flux above 0.1 GeV, which is the 
effective energy threshold for the EGRET blazar studies, and the effective cutoff seen at ~ 10 
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cm 2 s 1 corresponds to EGRET's quoted flux sensitivity. The solid line shows the expected flux 
distribution of blazar 7-ray sources under the assumption that the radio and 7-ray fluxes are 



linearly correlated [26, p7|]. The dashed vertical line is placed at the estimated flux sensitivity of 
AMS/7, which is about a factor of two lower than that of EGRET. Note that we have previously 
defined our flux sensitivity in terms of the integral flux above 1 GeV (not 0.1 GeV), since our 
detector's energy threshold is effectively at 0.3 GeV. To place the "equivalent" AMS/7 sensitivity 
on this figure with the EGRET data, we have assumed a nominal E~ 2 differential source energy 
spectrum. 

From the figure, it is seen that AMS/7 during its three-year mission should detect somewhere 
between 120 to 170 blazars, thus increasing the population statistics of these sources. This result 
is relatively model- independent; if we assume a spatially homogeneous population of sources 
in Euclidean space with no evolution, the count N versus flux S distribution follows a A(> 
S) oc 5~ 3 / 2 power law, as given by the dotted line in Figure 10, independent of the population's 
luminosity distribution function. This curve differs only slightly from the evolution-dependent 
curve at the AMS/7 sensitivity. 

With these increased statistics, improved 7-ray luminosity distribution functions can be de- 
termined e.g.,], and compared to those at other wavelengths, particularly at the radio [pS| , 
e.g.,]. Examination of the correlation between the radio and 7-ray luminosity functions would 
help, for example, to discriminate between 7-ray production models which predict different beam 
profiles of the emitted radiation (the 7-ray intensity as a function of viewing angle from the jet 
axis). For example, if the photons upscattered by the synchrotron-emitting electrons are "exter- 
nal" photons, coming from an accretion disk, the beam profile will be substantially narrower than 
if the upscattered photons are the synchrotron photons 



4.1.2 Detection of the Extragalactic Soft Photon Background 

All of the blazars detected by EGRET with high (> 6a) significance are found to have power law 
spectra, dN/dE cx E~ a , where the spectral indices a range from 1.7 to 2.5, with a mean index of 
2.15 ||. The measured energies range from the EGRET threshold of 0.1 GeV up to several GeV, 
the value of the upper energy being statistics limited and hence a function of source strength. 
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Under the assumption that one can extrapolate these power law spectra to TeV energies, several 
EGRET blazars should be detectable with ground-based atmospheric Cerenkov telescopes. In 



observations of EGRET blazars by the Whipple Observatory [31], however, no EGRET blazars 
other than the BL Lac object Mrk 421 were detected. For the blazars 3C 279, 1633+382, and 
2022-077, for example, the extrapolated EGRET flux is well above the sensitivity limit of the 
Whipple, yet these sources are not seen. Clearly something is happening to the 7-ray spectra of 
these blazars between a few GeV and the Whipple threshold energy (400 GeV at the time of these 
measurements) . 

Figure |ll| shows the EGRET power law spectra for four blazars extrapolated to 700 GeV. 
It should be noted that the majority of EGRET blazars are observed to be variable sources, 
with time- varying 7-ray flux amplitudes and spectral indices ||. We have somewhat arbitrary 
chosen a single flux amplitude and spectral index for each source from the Second Catalog p3] 



for purposes of illustration in Figure 11. Also shown are 7-ray flux upper limits for 1633+382 and 
3C279 measured by the Whipple ACT, which are seen to be far below the extrapolated power 
law spectra for these sources. 

There are a few possible reasons for a spectral cutoff between a few GeV and a few hundred 
GeV: (1) a cutoff in the electron spectrum due to synchrotron and inverse-Compton cooling may 
cause a corresponding cutoff in the upscattered 7-ray spectrum (an "intrinsic" cutoff mechanism) 



[33, e.g.,]; (2) 7-ray extinction through pair-producing interactions with the diffuse infrared (IR) 
to ultraviolet (UV) photons of the intergalactic photon background (IB) (an "extrinsic" cutoff 
mechanism) j34], ^2|| ; and/or (3) 7-ray extinction through pair-producing interactions with photons 
in the local environment of the jet, produced either within the jet or from the neighboring accretion 
disk (an "intrinsic" cutoff mechanism). 

By directly measuring the spectral cutoffs of the stronger blazar souces such as 1633+382 and 
3C279, AMS/7 may determine whether these known cutoffs are due to intrinsic mechanisms, or 
to the extrinsic process of interaction with IB photons. The main criterion for identifying the 
"extrinsic" mechanism is that the cutoff energies should decrease with increasing source redshift 
z, since the integrated opacity of the IB increases with source distance. The opacity r to a source 
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of redshift z e at a measured 7-ray energy of Eq is given by 

, , [ z ° dt f 2 x f°° . . .3^(2) . , , . 

T(E ,z e ) = c dz— dx- dv(l + z) — — o" 77 (s), (9) 
Jo az Jo 2 Jo hv 

where u v is the differential IB energy density at frequency v, c 77 (s) is the 77 — > e + e~ cross 

section as a function of the center-of-mass energy y/s, t is time, and x = (1 + cos#), 9 being 



the angle between the 7-ray and the soft photon momenta. Figure |11] shows the calculated 7-ray 
attenuation for two different models of the near-IR to UV intergalactic photon background [SjJ; 
one sees that the cutoff energy decreases with increasing redshift of the source. This figure also 
shows the "fauxdata" points of the AMS/7 spectra for these sources acquired over a two-year 
period of operation, under the assumption (for illustrative purposes) of power law spectra that 
are constant in time. These are obtained by convolving the source spectrum with the instrumental 
energy resolution function shown partially in Figure |9| The energy bin widths are indicated in 
the fauxdata for 3C273; the statistical errors are also shown at each energy point for the four 
sources. The effects of attentuation (according to the solid line predictions) have been included 
in the calculation of the fauxdata. 

Cutoffs are evident in the AMS/7 spectra of the stronger sources at higher redshifts (3C279 
and 1633+382), at energies ranging from 30 to 100 GeV. At higher energies the statistics are 
too low (and the AMS/7 energy resolution is significantly worsened) for AMS to measure blazar 
cutoffs. Since the threshold condition for pair production is 2Ee(l — cos#) > Arrive 4 , where e is 
the background photon energy, for a 7-ray of energy -E/reV only those background photons with 
energy e e v > 0.3/i?TeV can contribute to the opacity. Thus AMS/7 wm om y be able to detect 
7-ray attentuation due to the IB of optical and UV photons. Therefore the observation of such a 
redshift-dependent cutoff would constitute a direct, finite measurement of the intergalactic optical 
and UV photon background, which has not yet been made. Conversely, the absence of a redshift 
dependence for observed blazar cutoffs would definitively point to an intrinsic cutoff mechanism 
in the 7-ray production in jets. 

Finally, it has also been pointed out that if the IB spectrum u v {z) were known to high precision, 
then a measurement of blazar cutoffs, should they be extrinsic, would provide a measure of 
Hubble's constant Ho [36|, since the dt/dz term in Eq. ^ contains a factor of Hq 1 . Unfortunately, 
it is unlikely that such accurate spectral information will be available in the forseeable future in 
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the optical or UV bands. 

4.2 Gamma Ray Bursts 

Gamma ray bursts (GRBs) are likely to be the most energetic phenomena that occur in the 
universe. These point-source events flare into being with short, intense bursts of x-rays and 7- 
rays, then quickly fade away. GRBs are isotropically distributed in the sky, with no statistical 
evidence for repeating bursts occuring from a single source. First detected in the 1960's by the 
Vela satellites (whose function it was to search for nuclear weapons testing), it is only within the 
last two years that this class of events has been identified as being cosmological in origin, implying 
energy outputs on the order of 10 51 to 10 53 ergs released within a time scale of seconds. This 
rivals the mean luminosity of the entire universe, approximately 10 53 erg/s p7| ]. 

One of the two most popular current models for GRBs is the merger of binary neutron stars 
(with a maximum energy release of 5 x 10 53 ergs) or of a neutron star and black hole companion, 
an event which is estimated to occur on the order of once every 10 6 years per galaxy, in excellent 
agreement with the observed rate of GRBs (about one per day by the BATSE instrument of 
CGRO) |3lfl . The contending model is that GRBs are the results of "hyp er novae" , extremely 
energetic supernovae that result from the gravitational collapse of very massive stars. 

The nonthermal x-ray and 7-ray radiation observed from GRBs is most likely produced by 
ultra-relativistic particles (electrons and protons) that are accelerated in the shock waves gener- 
ated in the expanding fireball that follows the progenitor event. The strongest detections of GRBs 
are in the x-ray region of hundreds of keV, in which the BATSE, BeppoSAX, ROSAT, and ASCA 
orbiting instruments operate. Some of the very strongest bursts have also been detected in GeV 
7-rays by EGRET (four bursts with 7-rays above 0.1 GeV), one GRB yielding photon energies 
up to 26 GeV |?{|. Given the relatively low sensitivity of EGRET to GRBs, the detection rate by 
EGRET is consistent with all GRBs having spectra that extend to GeV energies (and possibly 
higher), possibly with nearly constant values of E 2 dN/dE [fiof . 

It is in fact the highest energy 7-rays that are the most constraining on the fireball models 
of GRBs; with the significantly larger sensitivity to > 10 GeV 7-rays that AMS/7 has compared 
to EGRET, the timing and spectral data of GRB 7-rays detected by AMS/7 could be invaluable 
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in elucidating the nature of the GRB fireball. Such data could also be applied to the question 
of whether the ultra- high energy cosmic rays might originate in the fireball shocks of GRBs pi]]. 
The fundamental issue involves the "compactness" of the GRB fireball, which is a dimensionless 
ratio of source luminosity to source radius. The higher the compactness, the greater the opacity 
for 7-ray emission, due to the large density of lower energy photons in the fireball that are above 
threshold for pair production. Fireball dimensions are causally limited by x-ray timing variations, 
whose time constants are sometimes as low as ~ 10 -3 s. Based on distance estimates from optical 
measurements, the luminosity of the fireballs can be inferred. These combine to form an average 
optical depth for 77 — > e + e~ of ~ 10 15 (-E/10 51 ergs)(5T/10 msec) -2 [38|. In spite of the enormous 
estimated opacities, 7-rays are observed from GRBs. This conflict can be resolved if the fireball 
is expanding with a very large Lorentz factor. Relativistic kinematics then both boost the 7-ray 
energy and observed luminosity of the source, plus reduce the time constants of timing variations, 
so that the apparent optical depth is approximately a factor r 4+2a larger than the actual value 
(a is the spectral index of the GRB spectrum). This argues for bulk Lorentz factors > 10 2 . The 
higher the observed 7-ray energies, the larger the fireball Lorentz factor must be to allow these 
7-rays to escape the fireball region 

In addition, there is a relationship between the redshift of a given GRB and the maximum 
energy 7-ray that can be received at Earth, due to the effect of extinction via pair production off 



the cosmic optical and UV photon background (see Sec. 4.1.2 ). No photons with energies above 20 
GeV, for example, have been detected from sources known to have redshifts above unity. Detection 
by AMS/7 of any 7-ray energies above ~ 20 GeV, coupled with an independent lower limit to the 
redshift of the GRB source via optical measurements, would place limits on the energy densities 
of the cosmic photon backgrounds if the redshift were on the order of unity or greater. 

The time distribution of the 7-ray events, relative to the prompt x-ray bursts that are observed, 
is also of great importance. In one GRB, GeV 7-rays were observed by EGRET for over one hour 
after the initial x-ray burst A number of scenarios have been proposed to explain this 

delayed emission, ranging from generation in more distant (external) shocks by the fireball, to the 
deflection of e + e~ pairs in the intergalactic magnetic field, followed by inverse Compton scattering 



that produces the 7-rays [43]. 
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It is not possible to predict with certainty how many GRBs will be detected by AMS/7. 
Given the comparable total exposures (> 10 11 cm 2 -s-sr for AMS) of EGRET and AMS/7 over 
the lifetimes of the two instruments, one could use the limited statistics of EGRET's four GeV 
GRBs as a naive estimate. In fact, it is likely that fewer would be detected by AMS/7 due to 
the higher effective energy threshold of 0.3 GeV for AMS/7, compared to 30 MeV for EGRET. A 
GRB similar to the one strong burst observed by EGRET, GRB940217 [j39|, would certainly be 
detected by AMS/7 if the burst were within a few tens of degrees from the orbital plane of the 
Space Station. Since the FOV of EGRET is smaller than that of AMS/7, an d since EGRET was 
not targeting GRB940217 when it occurred, it is safe to say that the a priori detection probability 
of a strong GeV burst by AMS/7 would be larger than for EGRET. 

4.3 The Extragalactic Gamma Ray Background at High Energies 

EGRET has recently measured the spectrum of the isotropic extragalactic 7-ray background 



(EGRB) pOfl, as shown in Figure 12. Originally detected by the SAS-2 satellite telescope at 



energies above 100 MeV [44], the origin of this background is still uncertain, and may be due to 
diffuse processes such as secondary production from cosmic ray collisions with intergalactic gas 
[45, 46], or, more likely, the accumulated emission of unresolved point sources, such as blazars or 
other classes of active galactic nuclei Kt], E8|. 



4.3.1 The Role of Blazars 

In particular, the hypothesis that unresolved 7-ray blazars are responsible for the EGRB has 
been examined, with conflicting conclusions (2^, pq ]. Stecker and Salamon (1996), by assuming 
a linear relation between the radio and 7-ray luminosity functions of blazars |27|, fit EGRET's 
blazar flux distribution (see Section ^4.1.1| ), and also predicted the unresolved blazar contribution 
to the EGRB, shown in Figure [12|. They were able to account for 100% of the observed EGRB, 
but predicted a sharp cutoff near 20-40 GeV due to the effects of 7-ray extinction off the IB 
optical and UV photons [32]. A more recent estimate of the 7-ray luminosity function of blazars, 



however, finds that unresolved blazars contribute only 25% of the total EGRB, thus indicating 
the need for additional (nonblazar) sources. Such a conclusion would not be inconsistent with an 
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EGRB whose spectrum continues as a power law to energies above 100 GeV. 

Whether blazars can or cannot contribute 100% of the EGRB may be resolved with the 
increased statistics at the highest energies that AMS/7 will obtain. If EGRET-type blazars do 
contribute 100% of the EGRB, then a cutoff musthe present somewhere below 100 GeV. Although 
the latest EGRET data shows no evidence of a cutoff in the EGRB spectrum, and is consistent 
with a pure power law (Eq.[|) EGRET's systematic errors make the current data somewhat 
inconclusive. AMS/7 m two years of operation should detect ~ 10 5 EGRB 7-rays in the interval 
0.5-1.0 GeV, and ~ 10 3 between 50-100 GeV; assuming that systematics can be kept to the 
several percent level, a spectral cutoff occurring within the 20-40 GeV region would be detectable 
by AMS/7. 

4.3.2 Instrumental Backgrounds to the EGRB Measurement 

Any measurement of the EGRB requires that instrumental background events be kept to a lower 
rate than the EGRB flux. To estimate the effects of background due to collisions of cosmic ray 
electrons, positrons, and protons with the AMS instrument, we divided the 7-ray spectrum into 
four bins per decade of energy (from 0.5 to 100 GeV) and required that we investigate backgrounds 
down to a level of 20% of the EGRB rate in each energy bin. For example, in the 25-40 GeV bin 
it will take AMS/7 2.0 x 10 5 seconds to obtain 5 EGRB 7-rays. In our Monte Carlo analysis, 
therefore, we threw 2.0 x 10 5 seconds' worth of cosmic ray electron, positron, and proton flux at the 
instrument to generate a false 7-ray background. These cosmic rays were thrown isotropically over 
a zenith angle range of 0° < 9 < 110°, where the largest zenith angle corresponds to the location 
of the Earth limb at an orbital altitude of 400 km. The energies were distributed according to 
known electron, positron, and proton energy spectra as discussed below. 

For cosmic ray electrons above 10 GeV, we used the compiled data presented in Muller and 
Tang (1987) [^9| to estimate the electron spectrum, and for energies below 10 GeV we used the 
data of Golden et al. (1994) ]50|] . We note that this latter set of data was obtained by a balloon 
flight in northeastern Canada, where the geomagnetic rigidity cutoff is well below 1 GV; since 
AMS during most of its orbit will be in regions with much higher rigidity cutoff (up to 15 GV 
near the equator), the use of the Golden etal data for the low energy electron spectrum is quite 
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conservative. Approximately 10 electrons in all were thrown at the AMS instrument, of which 
0.1% met the basic reconstruction criteria for a 7-ray candidate event. Additional quality cuts 
were applied so that this background was reduced to at most one count per energy bin (for a total 
of 2 background events). These cuts were then applied to the actual 7-ray data as described in 
Section EO . The positron flux was assumed to be 0.1 times that of the electron flux over the full 
spectrum (measurements give this ratio as varying between 0.05 to ~ 0.2 above 0.5 GeV ^T|). For 
a given 7-ray energy bin [E\ ow , -EhighL onr y electrons and positrons of energies greater than E\ ow 
were considered to be potential sources of 7-ray background. This assumption, justified by the 
fact that none of the reconstructed background events had energies higher than their progenitor 
electrons or positrons, reduced computational demands significantly. 

For cosmic ray protons we used the recent absolute flux measurements of Menn etal (1997) [^] 
which were obtained in a balloon flight at low geomagnetic cutoff. As with the electrons, therefore, 
our estimated low energy proton flux is conservatively high. Unlike electrons, which can generate 
bremsstrahlung photons with nearly the full primary electron energy, protons collisions generate 
7-rays primarily through 7r° decays (direct 7-ray production is only at the level of 10% that of 
7r°s at 500 GeV/c, for example |53|]), where the secondary ir°s carry only a fraction of the primary 
proton energy. Therefore, even though the CR p/e ratio is ~ 10 2 , we do not need to follow 10 9 
protons (compared to 10 7 electrons) through the detector. Instead, when estimating background 
production for the 7-ray energy bin [-Ei OW ) £high]> only protons of energy > r]E\ ow , r] > 1, were 
considered. We chose the value 77 = 7 based on calculations of galactic 7-ray production via CR 



proton interactions with the interstellar medium by Mori [54]. Using the HADRIN, FRITIOF, 
and PYTHIA particle collision codes, and convolving with the CR proton spectrum, Mori found 
that for a fixed 7-ray energy EL roughly half the photons were produced by protons of energy 
> 7Ej. This reduced our computational load by ~ r/~ 17 given the integral proton spectrum of 
E~ 17 above a few GeV; approximately 3 x 10 7 protons were required for these studies. We found 
that with the quality cuts introduced to eliminate the electron-induced background events, we had 
at most one proton-induced background 7-ray per source-energy bin (for a total of 5 background 
events); all these had reconstructed energies of < 1 GeV. 
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4.4 Detection of Neutralino Annihilation 



It is well known that the luminous matter within galaxies, and within our Galaxy in particular, is 
far below that required to explain the gravitational dynamics of their components. The inability 
of the estimated mass from disk and gas components to account for the observed rotation curves 



of spiral galaxies, for example, is exemplified in Figure 13, implying the presence of a "dark 
matter" component whose density falls off roughly as r~ 2 , where r is the distance from the 
center of the galaxy. Although MACHOs (massive compact halo objects) may possibly account 
for a significant fraction (20% to 100%) of the inferred dark matter in our Galaxy |55| . |56[| , the 
role of MACHOs in the mass budget of the Galaxy is still very poorly defined, and may prove 



to be minor [57]. An alternative is that this dark component is comprised of WIMPs (weakly 
interacting massive particles), one candidate for which is the supersymmetric neutralino, x- (The 
subject of supersymmetric dark matter has recently been comprehensively reviewed by Jungman, 
Kamionkowski, and Griest (1996) J58| .) 

Should supersymmetric neutralinos be the dominant dark matter component, the mutual an- 
nihilation of these Majorana particles near the center of the galaxy would produce a number of 
signatures, one of which would be the emission of both line and continuum 7-rays whose charac- 
teristic spectra would distinguish them from other contributions to the diffuse 7-ray background. 
The flux of monochromatic (line) 7-rays per solid angle from XX ~> 77 or XX ~^ annihilations 
in the galactic halo is given by [^] 

dF {a xx v) 



dVt 



8vrM2 J 



P A Ji,ip)dm, 



(10) 



where p x is the neutralino mass density, M x is the neutralino mass, a xx is the annihilation cross 
section into 77 or Z°j, the leading factor r/ 7 = 2 for XX ~~ 77 an d % = 1 for XX ~^ 
is the angle between the viewing direction and the direction to the center of the Galaxy, and I 
is the distance along the viewing direction. Because the annihilating particles are identical, the 
production rate per unit volume is rj~ f (n x /2) (a xx v), where the factor of 1/2 in eliminates 
double-counting of annihilating pairs. Equation III] can be expressed in the form 
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where 

■w-f S'GbK (12) 

Estimates of the distance of the Sun to the Galactic center, Rq, vary from 6 to 10 kpc, with a 
weighted mean of 8.0 ± 0.5 |6(| (to be compared with the "standard" value of 8.5 kpc adopted by 
the IAU in 1985 |6l|| ). We take i?o = 8.0 kpc. The local dark matter density po is obtained from 
measurements of the circular rotation velocity v c (R) of our Galaxy, where R is the galactocentric 
radius. The contributions to v c {R) from the disk and bulge (baryonic) components are subtracted, 
leaving a residual of v c ^ a i {R) 2 ~ GM(R)/ R. Fitting this to various dark matter halo profiles 
then yields the parameter pq. This has been done recently by Dehnen and Binney p2fl , who obtain 
values of po ranging from 0.1 to 0.7 GeV/cm 3 , depending on the assumed halo density profile. 
These authors also conclude that the mass distribution of the Galaxy is currently ill-determined, 
so no one model is clearly preferable over others. 

The form of the halo density profile p{R) is critical to the question of observability of 7-rays 
from neutralino annihilation. The simplest assumption, that of an isothermal halo, leads to a 
density profile p which is constant for R <^ a, where o is the "core" radius, and which falls as r~ 2 
for R 3> a |63[ ], accounting for the flat rotation curves of spirals at large radii. High resolution, 
N-body numerical simulations of the collapse of cold dark matter (CDM) halos, however, show no 
signs of the existence of a halo "core" down to the smallest resolvable scales. Studies by Navarro, 
Frenk, and White (NFW) (1996,1997) @, corroborated by others f6§, show that CDM halos 
over a wide range of total mass obey a "universal" density profile, given by 

P{R) * R(l + R/ar (13) 
where a is a scale radius. Although this profile results in a singular annihilation rate at the galactic 
center, the total 7-ray flux over any finite solid angle about tp = is finite. Even steeper radial 
dependences have been proposed: Berezinsky, Gurevich, and Zybin (1992) ]67|] argue theoretically 
for p(R) oc R 8 down to a cutoff of R ~ 1 pc, but this result has been criticized as arising from 
their use of unphysical initial conditions [S8|. However, in the mass modelling of the Galaxy by 
Dehnen and Binney JS^ ], one of their models has an R~ 1S density profile down to a core radius of 
a = 1 kpc, and fits the dynamical data with a good x 2 - Recent analyses of dwarf and low-surface 
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brightness galaxies, whose dark matter fraction may reach levels of 95%, indicate the presence of 



central cusps with density profiles p(R) oc 1/i? 7 , 7 ~ 0.2 — 0.4 [69], significantly shallower than 
that of NFW. It is not clear, however, that these shallower cusps are present in larger galaxies, 
such as our Galaxy, since the cusp steepness may depend on how dynamically "relaxed" the halo 
has become ]69[ ]. 

To estimate the potential sensitivity of AMS to monochromatic photons from neutralino an- 
nihilation, we use three of the halo profiles determined by Dehnen and Binney, corresponding to 
their models 2d, 4d, and 2f, plus a generic quasi-isothermal profile. All the Dehnen and Binney 
halo profiles take the form 

where g = (R 2 + z 2 /q 2 ) 1 / 2 , R being the radial distance along the plane of the galaxy, z being the 
height above the plane, and q = 0.8 being the assumed oblateness of the halo. The halo profile 
exponents 7^ and j3h vary with the mass model, and are given in table below: 



model 


lh 


fa 


a (kpc) 


Pa (M / P C 3 ) 


po (GeV/cm 3 ) 


D&B 2d 


1 


3 


21.8 


0.006159 


0.34 


D&B 4d 


1 


3 


5.236 


0.1101 


0.43 


D&B 2f 


1.8 


1.888 


1.0 


0.4179 


0.31 


isotherm 






1.0 




0.30 



In models 2d and 4d, the density profile is fixed to follow the "universal" profile found by 
NFW, for which 7^ = 1 and (3 n = 3. The last halo profile is an approximation to the true 
isothermal profile (see Binney & Tremaine, 1987): 

'a 2 + R 2 \ 



p(R) = po 



+ R 2 



(15) 



Figure 14 shows the resulting J(V') curves for the various halo profiles. 

The remaining two parameters that determine the 7-ray flux are, from Eq. |ll|, the neutralino 
mass M x and the annihilation cross sections a xx that result in 7-ray production. Monochromatic 
photons from XX ~~ > 77 an d XX ~^ ^°7 are necessarily produced via box diagrams 71 1, and so 
have a smaller branching ratio that tree-level annihilation into W and Z° pairs. Secondary 7-rays 
from the decay of tree-level W's and Z 0, s give a much larger continuum flux, which, depending on 
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the SUSY parameters, may be well above the extragalactic 7-ray background flux (Section |4.3| ) 
near the Galactic center [74]. 



4.4.1 AMS Sensitivity to Monochromatic Gamma Rays from Neutralino Annihila- 
tion 

Monochromatic 7-ray detection is limited both by detector energy resolution (Fig. ||) and aper- 
ture. Should the expected event rate be significantly greater than unity, the line width must be 
small enough so that the signal is not lost in the galactic and extragalactic 7-ray backgrounds. 
The rapid deterioration of energy resolution above 100 GeV limits AMS detection of monochro- 
matic 7-rays to neutralino masses M x < 200 GeV. (This limit does not apply to detection of 
continuum photons from W and Z° decays; see Section 4.4.2| .) However, recent analysis by the 



DAMA/Nal Collaboration has confirmed the presence of an annual modulation in their detector's 
event rate (at 98.5% C.L.), suggesting the existence of a local WIMP density with a la mass range 
M x = 59±u GeV > and WIMP-nucleon scalar elastic cross section £a£™j^) = 7.0±?;| x 10~ 9 nb, 
where £ is the fraction of the local dark matter density that is due to WIMP neutralinos ]ff2]| . 
Within this mass range our rms energy resolution is flat at 24%. These mass and cross section 
ranges have also been shown to be fully compatible with a MSSM neutralino being the major 



component of dark matter in the universe [73]. 



Figure 15 illustrates the sensitivity of AMS to detection of line photons from the annihilation 
channel XX -* 77 as a function of halo model. Plotted are integrated 7-ray counts N(< ifj) 
as a function of angle ip from the Galactic center for the four halo profile models, accumulated 
over a two-year period by AMS, assuming M x = 50 GeV and a velocity-cross section product 
vo~ = 10~ 29 cm 3 /s. With an energy acceptance width of 20 GeV, approximately 80% of the line 
photons are captured. The dot-dashed line shows the galactic and extragalactic background count 
for this energy window; the parameterization of Bergstrom, Ullio, and Buckley was used for 
the galactic background contribution. The dashed and dotted lines correspond respectively to the 
"Berezinsky-like" halo (model 2f of Dehnen and Binney p2|) and the quasi- isothermal halo, while 
the upper and lower solid lines correspond to models 4d and 2d of Dehnen and Binney, which are 
the NFW-like models. 
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Bergstrom and Ullio [71] have recently performed a full one- loop calculation of the XX ~^ 77 



annihilation reaction, assuming the minimal super symmetric extension of the standard model 
(MSSM). They performed extensive scans in SUSY parameter space to map out the range of 
values of the annihilation cross section versus neutralino mass. We have reproduced their Figure 



3a in Figure 16, along with the AMS sensitivity for the two most optimistic halo profiles, models 
2f and 4d, expressed as a lower limit on (va) versus M x for detection of neutralino annihilation. 
The detection criteria were arbitrarily chosen to be iV 7 > 10 and iV 7 / ^JB~ > 6, where N~ is the 
total number of annihilation 7-rays detected over a 2-year operating period within the central few 
degrees of the Galactic Center, and £L is the total number of background photons from diffuse 
galactic emission detected within the energy window and the same solid angle about the Galactic 
Center. 

Ullio and Bergstrom fFOfl also performed a MSSM full one-loop calculation of neutralino an- 
nihilation via XX ~* 7^1 assuming a pure higgsino state for the neutralino. For this channel the 
7-ray has an energy Ey = M x — m 2 z /4:M x . In Figure |l7] we have reproduced Figure 6 of Ullio 
and Bergstrom, which shows their vaz-y versus higgsino mass (solid line), along with the unitary 
lower bound on this cross section (dashed line). Superimposed on this plot is a line segment 
which shows the minimum vaz-y required for detectability by AMS over a 2-year period, assuming 
the "Berezinsky-like" halo model 2f of Dehnen and Binney [p^] . None of the other halo models 
intersect the Ullio and Bergstrom curves. Thus for this pure higgsino annihilation channel only 
the most singular halo profile can yield a detectable signal. 

In general, then, only the more singular of the halo profile models, in which the dark matter 
density continues to increase as i2~ 7 , 7 ~ 1 — 2, can give hope for line photon detection by AMS, 
and then only in a restricted region of MSSM parameter space. It should be kept in mind that a 
fraction of this restricted region may already have been eliminated by the absence of a detection 
by EGRET, but such an analysis has yet to be completed |fT5|] . 

4.4.2 AMS Sensitivity to Continuum Gamma Rays from Neutralino Annihilation 

Continuum photons from the decay of Ws and Zs produced by tree-level annihilations XX ~> 
WW, ZZ are far more abundant than the line 7-rays considered above, and so possibly provide 
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the greatest chance for detection, assuming the neutralino mass is above that of the W and Z. 



Berstrdm, Ullio, and Buckley [74] have used the PYTHIA Monte Carlo code [76| to determine 

the the continuum photon energy distribution from XX annihilation into WW, ZZ, and qq pairs. 

They find that the hardest distributions come from WW and ZZ pairs, and that the annihilation 

rate is dominated by WW and ZZ production at higher (higgsino-like) neutralino masses. From 

a scan of MSSM parameter space, they find that the models with the highest rate of WW and 

ZZ production yield an observed differential photon flux of 

dF y in / 300 GeV\^ dNy . . a i i« , r i 
— 1 ~ 1.2 x 10~ 10 — — TEr J W cm" 2 s _1 sr- 1 GeV _1 , (16) 



dE 7 \ M x J dE^ 



where 



dx x dE 1 x 1 - 5 ' v ; 

for M x > 300 GeV. The J(ip) factor is given by Eq. 0. 

Using the J(ip) functions of the four halo models described above (Fig. 14), the differential 



flux coming from a solid angle Q about the Galactic Center can be obtained and compared to 
the flux from the galactic and extragalactic backgrounds. Using the collection area*time product 



At(E^,5) (Fig. |6|) for AMS, the differential event count is obtained. Figures 18 to 20 show these 
event counts after integration over ip out to three different solid angles about the Galactic Center, 
ft = 10~ 3 , 10 -1 ' 5 , and 1.0 steradians. The dot-dashed line is the diffuse galactic and extragalactic 
background photon count; the upper and lower solid lines correspond to halos for models 4d 
and 2d, respectively; the dashed line corresponds to the "Berezinsky-like" halo model 2f; and 
the dotted line comes from the quasi-isothermal halo model. The faux-data points shown in the 
figures take into account the effect of detector energy resolution; these are the counts obtained in 
each energy bin after convolution with the Monte Carlo energy resolution functions (Fig. ||). 

There are a number of salient features, the primary one being the strength of the signal 
compared to the diffuse background for essentially all of the halo models considered. For this 
particular (optimized) MSSM model it is likely that EGRET would have already provided a 
detection given any of the halo models considered here. Although there is a large variation in 
signal between the various halo models, it appears nevertheless that for any of these halo models 
a substantial region of MSSM parameter space can be constrained by the absence of a continuum 
7-ray signal. 
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These figures also show that the differing angular dependences of the background and signal 
result in there being an optimal solid angle f2 that maximizes signal to noise, with more sin- 
gular profiles having a stronger signal to noise at smaller solid angles. These differing angular 
dependences can also be used to discriminate the signal from the diffuse backgrounds. 

Finally, we point out that continuum photons result from XX ~^ bb as well as from WW and 
ZZ pair production, producing scaled spectra similar to those from W and Z decays. This implies 
that AMS/7 sensitivity to continuum 7-rays may extend down to the current LEP limits on the 
neutralino mass jr7j . 

5 Conclusions 

We have shown that with minor modifications the Alpha Magnetic Spectrometer can become a 
powerful 7-ray detector as well, with overall performance characteristics being comparable, if not 
superior, to those of EGRET. With 7-ray energy resolution extending past 100 GeV, and with an 
aperture that is nearly flat above ~ 3 GeV, AMS/7 can address a number of outstanding issues in 
7-ray astrophysics that relate to the relatively unexplored region of Ey = 20 — 200 GeV. For one, 
AMS/7 wm likely confirm or refute the hypothesis that unresolved blazars are responsible for the 
bulk of the extragalactic 7-ray background; AMS/7 wu l a ^ so extend the spectrum of the diffuse 
galactic background to above 100 GeV, helping to resolve current difficulties in interpreting the 
EGRET diffuse galactic background measurement ]7q| . 

AMS/7 should roughly double the total number of blazars detected in 7-rays, and will be 
enable multiwavelength observational campaigns to include the GeV region of blazar spectra 
during the flight years of 2003-2006. There is an additional possibility that an indirect detection 
of the cosmic UV and optical photon background can be made through the detection of extinctions 
in high-redshift blazars above ~ 20 GeV. AMS/7 wm a lso likely observe GeV 7-ray emission from 
one or more gamma ray bursts during its operational lifetime. 

AMS/7 wu l a l so search for both line and continuum emission of 7-rays from the region of the 
Galactic Center created by the annihilation of dark matter neutralinos. Although the sensitivity to 
line emission appears marginal, there is nevertheless a finite, though small, region of halo/MSSM 
phase space which allows a detection by AMS. However, a much larger region of dark-matter 
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halo/MSSM parameter space can be constrained in a search for continuum 7-rays by AMS/7. 

Finally, we note that even higher sensitivities can be reached with the addition of a high 
granularity calorimeter below the magnet (option (b)); this will be the subject of future work. 
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Figure 1: A cross section of the AMS instrument. 
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Figure 2: The momentum resolution of the AMS spectrometer for e, 7T, K, p, 3 He, 4 He, and 12 C. 
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Figure 3: Angular resolution of AMS/7 as a function of primary 7-ray energy (filled squares), and that 
of GLAST H) (dashed line) in the energy interval 0.3 to 100 GeV. 
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Figure 4: (a) Histogram of detected 7-rays as a function of incident zenith angle 9 for 10 GeV 7-rays. (b) 
Effective 7-ray detection area A(E, 9) versus zenith angle at 10 GeV. 
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Figure 5: AMS/7 aperture as a function of 7-ray energy. The statistical errors are comparable in size to 
the circular data points. 
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Figure 6: Accumulated area*time product over one year of AMS operation as a function of declination; 
each curve is marked by the 7-ray energy (in GcV). The declination of the Galactic Center (GC) is also 
indicated as a dashed vertical line. 
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Figure 7: A 3-D plot and 2-D contour plot of AMS/7's point source sensitivity n versus celestial coor- 
dinates. Recall that no is the minimum amplitude of the source's differential flux at 1 GeV required for a 
5a significance detection. The units for n in the 3-D plot are 10~ 8 cm _2 -s _1 -GeV _1 . Large photon fluxes 
from the diffuse galactic background are responsible for the deterioration of sensitivity near the Galactic 
plane. 
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Figure 8: The reconstructed energy distributions for four different primary 7-ray energies, 0.3, 3, 20, and 
100 GeV. 
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Figure 9: The mean fractional energy loss (open square), rms energy resolution (filled square), and width 
of the "gaussian" peak of the distribution (filled circle) as a function of primary 7-ray energy. 
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observed integral flux (>100 MeV) (10 8 cm" 2 -s" 1 ) 

Figure 10: Full sky blazar count as a function of the detected integral 7-ray flux above 0.1 GeV. The 
ordinate gives the total source counts within one-fifth-decade intervals of integral flux. The unfilled circles 
represent EGRET detections, taken from their Second Catalog [^|; the solid line is the predicted sky count 
from Stecker and Salamon (1996), and the dotted line is the Euclidean relation, iV(> S) oc S~ 3 ^ 2 . The 
dashed line shows the point source integral flux sensitivity of AMS/7. 
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Figure 11: Power law spectra of the blazars 1633+382, 3C273, 3C279, and Mrk 421, extrapolated from 
the measurements of EGRET. The redshifts z of these sources are indicated in the figure. The dotted lines 
show the pure power law extrapolation, while the solid and dashed lines show the effect of 7-ray extinction 
via pair production off the intergalactic soft photon (IR to UV) background according to two different 
models of the soft photon background [j32[ . The open triangle and circle represent 7-ray flux upper limits 
from 3C279 and 1633+382, respectively, as measured by the Whipple Observatory These are seen to 
be well below the naive extrapolation of the EGRET power law spectra of these two sources. Note that we 
show only one power-law representation for each source; some sources, such as 3C279, are highly variable, 
with significant shifts occurring in both the integrated 7-ray flux and spectral index during flare states 
Pi- Also shown are the expected AMS/7 spectrum measurements corresponding to a two-year integration: 
3C279 (filled upper triangles), 1633+382 (filled circles), Mrk 421 (filled lower triangles), and 3C273 (filled 
squares). 
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Figure 12: The extragalactic gamma ray background as measured by EGRET (filled circles), along with 
their fitted (straight solid line) power law spectrum. Also shown are two versions of a theoretical spectrum 
calculated under the assumption that 100% of the EGRB is due to unresolved blazars 26| 32|, which includes 
the effects of 7-ray attenuation via interactions with the IB optical-to-UV photons. (The solid and dashed 
line curves differ in assumptions regarding the effects of evolving cosmic metallicity on the optical and UV 
photon backgrounds.) The peaks in E 2 dN/dE at ~ 20 GeV in these theoretical spectra are sensitive to 
assumptions regarding the spectral index distribution of blazars, which appear to be inconsistent with the 
EGRET data, but could be readily modified to be more so. The one feature which is not malleable is the 
spectral cutoff at ~ 20 — 30 GeV; for an EGRB of blazar-dominated origin a spectral cutoff in the 20-50 
GeV region must exist. Should the EGRB spectrum remain a power law out to 100 GeV, the blazar-origin 
hypothesis would be ruled out. The error bars on the EGRET data (filled circles) are mostly systematic 
(13% for 0.03 to 10 GeV data; 30% above 10 GeV) [Sreekumar etal 1998], representing over 3 years of 
EGRET data. AMS statistical errors, based on one year of operation, are smaller than the data point 
symbols (filled triangles). For example, between 50 to 100 GeV, AMS should see 1.3 x 10 3 EGRB 7-rays 
in one year, assuming the EGRET EGRB spectrum. 
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Figure 13: The rotation curve of the galaxy NGC 6503 is shown as a function of radius from the galaxy 
center, as measured by Doppler shifting of the 21-cm radio line from gaseous HI in the disk. The flatness 
of the rotation curve at large radii implies an integrated galaxy mass M(r) that falls off roughly as r _1 , 
inconsistent with the observed mass profiles of the visible disk and gas components of the galaxy. The 
inclusion of a dark matter halo component with density profile ~ r~ 2 provides an excellent fit to the 
observed rotation curve. This figure is taken from Begeman, Broeils, and Sanders (1991) [Q. 
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Figure 14: The integrated 7-ray flux J(ip) versus the angle ip from the Galactic center for the tabulated 
halo dark matter density profiles, where J(ip) is defined in Eq. [l^. The dashed, upper and lower solid, and 
dotted lines correspond respectively to D&B models 2f, 4d, 2d, and the quasi-isothermal halo profile. The 
most optimistic profile is model 2f, which is very close to that of Berezinsky et al (1992). Models 4d and 
2d are both NFW halo profiles. 



47 




-4 -3 -2 -1 

10 10 10 10 1 

\|/ (radians) 



Figure 15: Integrated 7-ray counts N(< tp) as a function of angle ip (in radians) from the galactic 
assuming M x = 50 GeV, producing a narrow line at E 7 = 50 GeV. See discussion in text. 



48 




20 30 40 50 60 70 80 90 100 

Neutralino mass (GeV) 

Figure 16: The velocity-cross section product required for AMS detection of a neutralino induced 7-ray 
line as a function of neutralino mass and halo profile model. The dashed line corresponds to D&B model 
2f, which is the "Berezinsky-like" model; the solid line corresponds to D&B model 4d, the more productive 
of the two "NFW-likc" halo profiles. These curves are superimposed on Fig. 3a of Bergstrom and Ullio 
(1997), which gives their calculated cross sections based on a scan of SUSY parameter space. 
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Figure 17: The velocity-cross section product required for AMS detection as a function of neutralino mass 
and halo profile model. These curves are from Fig. 6 of Ullio and Bergstrom (1997), which gives their 
calculated cross section for the single-photon production channel \X ~ * 1% or a pure higgsino state. The 
dashed-line line segment corresponds to D&B model 2f, which is the "Berezinsky-likc" model. The other 
halo profiles do not intersect the Ullio and Bergstrom curves. (The solid line segment in the upper-left 
corner corresponds to D&B model 4d.) 
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Figure 18: The differential continuum 7-ray event count from neutralino (higgsino) annihilation in the 
Galactic halo detected by AMS over a 2-year time period, within a solid angle of 10 -3 sr about the Galactic 
center, based on the cross sections for \X ~ > WW, ZZ of Bergstrom, Ullio, and Buckley Jrl| . The dashed, 
upper and lower solid, and dotted lines correspond to the various halo models: Models 2f, 4d, 2d, and 
quasi-isothermal, respectively. The dot-dashed line shows the sum of the integrated diffuse galactic and 
extragalactic 7-ray background 7-rays detected by within the same solid angle. The faux-data points 
show the energy bin counts obtained by AMS after convolution with the energy resolution function of the 
detector. (Faux data points are not shown for the quasi-isothermal halo model.) 
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Figure 20: The same as Figure ttq, except the fluxes are those that arrive within a solid angle of 1 sr 
about the Galatic center. 
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